Mon. Not. R. Astron. Soc. 000, ITHHl (2011) Printed 9 July 2012 (MN M£X style file v2.2) 



Models for the lithium abundances of multiple populations 
in globular clusters and the possible role of the Big Bang 
lithium 

O. Francesca D'Antona 1 , Annibale D'Ercole 2 , Roberta Carini 1 ' 3 , Enrico Vesperini 4 , 
^ & Paolo Ventura 2 * 

r ^ " 1 INAF, Osservatorio Astronomico di Roma, Via Frascati 33, 1-00040 Monteporzio Catone (Roma), Italy. 
t~ 5 , 2 INAF- Osservatorio Astronomico di Bologna, via Ranzani 1, 1-40127 BOLOGNA (Italy). 

3 Department of Physics, Universita di Roma "La Sapienza" , Roma, Italy 

4 Department of Physics, Drexel University, Philadelphia, PA 19104, USA 



Accepted . Received ; in original form 



o 

CU ABSTRACT 

6 

1 Globular cluster stars show chemical abundance patterns typical of hot-CNO pro- 

I cessing, and photometric evidence for the presence of multiple populations. Lithium 

C$ ■ is easily destroyed by proton capture in stellar environments, so determining its abun- 

dance may be crucial to discriminate among different models proposed to account for 
the origin of the gas from which the multiple populations form. In order to reproduce 
the observed O-Na anticorrelation and other patterns typical of multiple populations, 
■s^j- " the formation of second generation stars must occur from the nuclearly processed stel- 

■^j- , lar ejecta, responsible of the chemical anomalies, diluted with pristine gas having the 

l/^ ■ composition of first generation stars. This gas is either a remnant of the first phases of 

star formation, or it has been re-accreted on the cluster core after the end of the SN II 
explosions. As a consequence, the lithium abundance in the unprocessed gas — which 
■ is very likely to be equal to the lithium abundance emerging from the Big Bang — 

\ affects the lithium chemical patterns among the cluster stars. This paper focuses on a 

scenario in which processed gas is provided by asymptotic giant branch (AGB) stars. 
We examine the predictions of this scenario for the lithium abundances of multiple 
populations. We study the role of the non-negligible lithium abundance in the ejecta 
of massive AGB (A(Li)^2), and, at the same time, we explore how our models can 
' constrain the extremely large — and very model dependent — lithium yields predicted 

by recent super-AGB models. We show that the super-AGB yields may be tested by 
examining the lithium abundances in a large set of blue main sequence stars in w Cen 
and/or NGC 2808. In addition, we examine the different model results obtained by 
assuming for the pristine gas either the Big Bang abundance predicted by the stan- 
dard models (A(Li)=2.6-2.7), or the abundance detected at the surface of population 
II stars (A(Li)=2.2-2.3), and we show that, once a chemical model is well constrained 
by a comparison with the observations, the O-Li distribution could perhaps be used 
to shed light on the primordial lithium abundance. 

Key words: stars: Population II; stars: abundances; stars: AGB and post-AGB; 
globular clusters: general. 
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1 INTRODUCTION 

The spread of light elements abundances in globular clus- 
ters is larger than in field stars of similar metallicity and 
requires the presence of a population of stars formed from 
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matter processed through the hot CNO cycle and by other 
proton-capture reactions on light nuclei (hereafter we will 
refer to this population as second generation, SG). The site 
of processing of this matter in the first generation (FG) 
stars of the proto-cluster is a subject of intensive investi- 
gation. The two main scenarios are: the "AGB scenario", 
in which the site of processing are the hot bottom burning 
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envelopes of massive asymptotic giant branch (AGB) stars 
dVentura et alj200ll ;|P' Antona fc Calo i 2004; Karakas et al.l 
2006), and the "fast rotating massiv e stars (FRMS) sce- 
nario" (IPrantzos fc Charbonnel 120061 ; iMevnet et all l200d : 
iDecressin et al. 2007al ). in which the site is the interior of 
fast rotating massive stars. 

Problems are present in both scenarios (see, e.g., 
iRenzinil [2008, for a critical review). In all the scenarios 
proposed so far, the composition of SG stars is, in most 
cases, best explained by dilution of the pollu ting ejecta 
with pristine gas (see, e.g. iD'Ercole et aljfeoill ). One pos- 
sible discriminant between the two scenarios is the com- 
parison between the abundance of lithium in FG and SG 
stars. In fact, massive stars destroy lithium, while massive 
AGB produce it at the beginning of the HBB (|Ventura et al.l 
2002), so possibly the prediction of the two models differ. 
Also for other scenarios lithium can be a powerful discrimi- 
nant: the polluting matter is Li-free also if it comes from 
runw ay collision between massive stars (jSills fc Glebbeekl 
|2010|). or non-conservative evolution of massive binaries 



( de MinketHI 120091 ). while the diluting gas may be Li 
free, if it comes from ma ss loss from FG stars by winds 
IIGratton fc Carretta|[2O10l ). and also if it is made up by the 
matter in non conservative evo lution of interacting close bi- 
naries l|Vanbeveren et al.ll2012|') or stripped dur ing close en- 
counters in the cluster core " ifCarini et al.ll2012l ). 

Self-consistent models to study the chemical evolu- 
tion of lithium in multiple-population GCs are still lacking, 
while some sets of obse rvational data are already available. 
IDecressin et al.l |2007bl ) sho w that the anticorrel ation Li- 
Na in the data of NGC 6752 jPasauini et al.ll2005l) are fully 
compa tible with a simple dilution model, but Ishen et all 
l|2010l ) present new data Li-O, probably not compatible 
with diluti on with Li-free matte r. In other clusters, like 
NGC 6397 ilLind et al.ll2009l ). M 4 ilD'Orazi fc Marindl2010l ; 
iMonaco et al.ll2012l ). 47 Tuc i|D'Orazi et al.ll2010l ) the FG 
and SG stars have very similar lithium content, and Li- 
depleted stars may be attributed to convective dilution. 
Naively, the similar abundances in these clusters seems in 
better agreement with the AGB scenario. 

In this paper, we study the lithium abundance patterns 
resulting from t he ch emical evolutio n models presented by 
ID'Ercole et~ai1 (|2010l ) (Paper I) and ID'Ercole et~ai1 (|2012T ) 
(Paper II). The models are based on the AGB scenario, 
and on the hydrodynamical computations by ID'Ercole et aU 
(2008), and can reproduce the chemical patterns both 
in clusters having an extended O-Na anticorrelation, like 
NGC 2808, and in clusters showing a mild anticorrelation, 
like M 4. 

The outline of the paper is the following. In Sect. [2] we 
briefly summarize the model, describe the lithium yields of 
the AGB and super-AGB models. As there is a still debated 
discrepancy between the predictions for primordial lithium 
from the standard Big Bang nucleosynthesis and the abun- 
dance observed in the atmospheres of population II stars, an 
additional hypothesis must be made, concerning the choice 
of the initial lithium abundance in the diluting p r istine gas. 

Sect. [3] compares the data bv lMonaco et al.l (|2012l ) for 
M 4 with the models that account for the O-Na anticorre- 
lation. 

In Sect.|4]we show that the predictions for lithium in the 
model for NGC 2808 is relevant also for a comparison with 



the data bv lShen et all l|2010l ) for NGC 6752. We will show 
that in this case the lithium yield of AGB stars is necessary 
to reproduce the observations, and a better determination 
of the O-Li patterns may provide an independent constraint 
on the lithium abundance emerging from Big Bang. 

Sect.[S]shortly examines the case of dilution of nuclearly 
processed ejecta with lithium-free gas. 

The possible constraints posed by multiple population 
abundance patterns on the Big-Bang lithium abundance are 
further discussed in Sect. [5] In Sect. [7] we summarize our 
conclusions. 



2 INPUTS OF THE MODEL 

2.1 The chemical evolution model 

The chemical evolution model, described in detail in Paper 
I and Paper II, assumes that the proto-GC becomes com- 
pletely devoid of pristine gas at the end of the SN II epoch, 
and begins a new star formation epoch when the low veloc- 
ity winds of the super-AGB stars first, and of the massive 
AGB stars later on, collect in a cooling flow in the cluster 
core. Formation may occur directly from the super-AGB 
and AGB ejecta, but a phase of re-accretion of pristine gas 
(see e.g. D'Ercole et al. 2008, 2011) may lead to dilution of 
the ejecta with the original matter. This phase is character- 
ized by an epoch of maximum accretion of pristine gas, t acc , 
and a timescale of accretion r pr i, that in our formulation is 
the standard deviation of the gaussian describing the pro- 
cess. A third parameter, the end (t on d) of the second star for- 
mation epoch, provides a further independent timing of the 
process. It is also possible to accumulate and mix the ejecta 
for a time tt before starting star formation. Further details 
on the general behavior of models and on the adopted sym- 
bolism can be found in Paper I. In Paper II we introduced in 
the model the new yields for super-AGB models from 6.5 to 
8Mq computed bv lVentura fc D'Antonal j201lf>. that com- 
plement the results bv l Ventura fc D'Antonal (J2009J) , to revise 
the results presented in Paper I. In this work we explore the 
multiple generation lithium distribution resulting from some 
of the models computed in Paper II for the clusters M 4 and 
NGC 2808. The main parameters of the chosen models are 
summarized in Table [I] We then extend the model discus- 
sion to other clusters for which data of lithium abundances 
are available (NGC 6397, NGC 6752, oj Cen). 



2.2 Lithium abundances in AGB and super-AGB 
ejecta 

For the sake of completeness, we list in Table H] the lithium 
abundances in the ejecta of AGB and super-AGB stars 
adopted in this work, and the evolutionary times of the 
corresponding masses. Notice that the assumptions on the 
initial lithium abundance in the gas forming the FG stars 
have no consequences on the yields. In fact, the initial 
lithium is fully depleted in the stellar envelopes, by the time 
the lithium produc tion due to the Cam eron-Fowler (1971) 
mechanism begins l|Ventura et al.ll2002l ) .There are some dif- 
ferences between the s e abu ndances and the ones listed in 
IVentura fc D'Antonal (|2010T ). due to recomputation of some 
evolutionary tracks. There is a discontinuity in the lithium 
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Table 1. Dynamical input parameters for the simulations 



Table 2. Lithium average abundances in the ejecta for Z=10 



Model 


^acc 


Tpri 


^end 


t{ 


P*,FG 


P0,pr 


V 


X 


M 4-0 


65 


2 


105 





9.4 


0.091 


0.1 


0.7 


M 4-2 


43 


2 


60 





940 


0.05 


0.1 


0.5 


M 4-3 


43 


2 


50 


45 


940 


0.04 


0.5 


0.5 


M 4-4 


43 


2 


55 


45 


940 


0.04 


0.5 


0.5 


M 4-5 


48 


10 


58 


48 


940 


0.08 


0.5 


0.5 


2808-1 


65 


8 


90 





240 


0.0095 


1 


0.4 



t acc : time at which maximum accretion of pristine matter occurs (Myr) 
r pr ;: timescale of the pristine gas accretion process (Myr) 
t cn( j: time at which the SG star formation ends (Myr) 
tf : time until which the AGB ejecta are accumulated before star 
formation begins 

p^FG 1 densiy of FG stars (Mqpc - 3 ) 

po.pr: densiy regulating the amount of available pristine gas, 
normalized to yg 
v: star formation efficiency 

x = P* SG/P*tot: ratio between the SG stars and the total 
nowadays alive stars (p* I tot = P*,FG + P*,Sg) 
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Figure 1. The full lines represents the lithium versus sodium 
(left), and lithium versus oxygen (right panel) abundances in the 
ejecta of super-AGB and AGB stars. Each dot represents an ini- 
tial mass, the maximum and minimum mass are labelled. The 
dash-dotted line represent the oxygen abundance that would be 
seen in giants having very large helium abundance (for stars di- 
rectly formed from the super-AGB ejecta), thanks to the 'in situ' 
deep mixing described in D'Ercole et al. (2012). Notice that we 
plot again the initial lithium abundances in the gas, but actually 
in these giants lithium is drastically reduced, both by standard 
convective mixing, and, even more, by the same convective deep 
mixing to which the oxygen reduction is attributed. 



M/Mq 


r/10 6 vr 




M„„„„ /Mri 
m core / lvA y 


core 


3.00 


332 


2.77 


0.76 


CO 


3.50 


229 


2.43 


0.80 


CO 


4.00 


169.5 


2.20 


0.83 


CO 


4.50 


130.3 


2.00 


0.86 


CO 


5.00 


103.8 


1.98 


0.89 


CO 


5.50 


85.1 


1.93 


0.94 


CO 


6.00 


71.2 


2.02 


1.00 


CO 


6.30 


65.2 


2.06 


1.03 


CO 


6.50 


61.5 


2.36 


1.09 


O-Nc 


7.00 


53.7 


2.12 


1.20 


O-Nc 


7.50 


46.8 


2.75 


1.27 


O-Nc 


8.00 


38.8 


4.39 


1.36 


O-Ne 



abundances as a function of the mass, that can be seen 
from Fig. [T] and Table [2] for increasing initial mass, the 
lithium abundances first decreases, reaching a minimum of 
A(Li)~1.9 at 5.5 Mq, then increases again. The initial de- 
crease is due to the increasing temperatures at the bottom 
of the convective envelopes, that imply a faster lithium pro- 
duction and destruction. For the larger masses, the higher 
luminosities provoke larger mass loss rate at the stage of 
lithium production, and the abundance increases again. Nev- 
ertheless, there is a second minimum in the abundance at a 
mass of 7 Mq, that is the first mass for which the lithium 
production occurs before the beginning of the thermal pulse 
phase. For M=7.5 and 8Mq the mass loss is dominant and 
huge quantities of lithium are lost to the interstellar medium, 
finally raising the lithium in the ejecta to values such as 
A(Li)=4.39 in the 8Mq. We warn that these numbers are 
dramatically d ependent on the adopted m ass loss rate (see, 
for a discussion I Ventura fc D'Antonal2010l ). nevertheless the 
discontinuity at 7 Mq can not be attributed to simple nu- 
merics, and comes out from the above physical explanation, 
so we keep it in the ejecta table. 

Fig. [T] summarizes the Li-Na-O abundances in the 
ejecta of super-AGB and AGB stars from Table 1 in Pa- 
per II. The most striking characteristic of these abun- 
dances is the very large Li and Na abundances reached in 
the ejecta of the 8Mq sta r. As remarked above (see also 
IVentura fc D'Antondl2010l ). this result is very uncertain, as 
it may be due to an overestimate of the mass loss rate during 
this evolution. Values of lithium abundances much smaller 
both for the 8 and 7.5 Mq stars are possible, without sig- 
nificantly affecting the sodium and oxygen abundances. As 
we will further discuss below, new observations will allow to 
put tighter constraints on the lithium yields of these stars. 

In Fig. [T] we also show the much smaller oxygen val- 
ues (due to deep extra-mixing) that Paper II assumes to 
be present at the surface of SG stars formed directly from 
the super-AGB ejecta, and thus having very high helium 
content. The values plotted refer to the abundances in the 
ejecta but, obviously, the lithium abundances in giants will 
be much smaller than those plotted since lithium will be 
strongly diluted, both by standard and deep mixing. The 
abundances in the ejecta of the masses 5-6.5 Mq that should 
contribute to the SG are A(Li)~2.0. Consequently, the AGB 
ejecta will be important for the final lithium content in the 
SG stars only in case B of Table [3] (see Sect. 12.31) . namely, 
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when the Li abundance in the diluting pristine gas is as low 
as that shown in the atmosphere of normal population II 
stars: if the initial Li is A(Li)=2.6, it is difficult that abun- 
dances 0.6 dex smaller (a factor 4 smaller!) in the ejecta 
may alter the Li content of the mixture in a significant way 
different from ejecta with no lithium at all. 



2.3 Assumptions for the Big Bang Lithium 
abundance 

The application of the dilution model for lithium abun- 
dances is not as straightforward as for the other ele- 
ments. In fact, we do not know observationally the lithium 
abundance of the pristine gas. The standard Big-Bang 
nucleosynthesis (BBN), with the cosmological parame- 
ter ?7=N(baryons)/N(photons) constrained by the obser- 
yations of the satell ite WMAP, provides A(Li)=2.64±0.04 
jSpergel et all l2007lh or even higher, A(Li)= 2.72±0.05 
when updated rates are taken in t o acc ount for the 
3 He(a,7) 7 Li reaction l|Cvburt et al.l 120081 ). The values 
computed from the analysi s of spectra of population II 
dwarfs are A(Li)=2.37±0.06 (|Melendez fc Ramfrezll2004h, or 
A(Li) =2.23±0.07 in the non-LTE analysis bv lAsplund et all 
(20060 We have to deal with the abundances into play by 
cons idering all th e possible explanations of this discrepancy 
(see lFieldsll201ll . for a recent review). The possible choices 
are: 

(i) w e reject the stand ard BBN scenario (for a review, 
see, e.g. Ilocco et al.ll2009h and assume that the abundances 
observed in pop. II stars (and in the FG stars of GCs) are 
the primordial abundances. Thus we adopt A(Li)=2.20 or 
2.30 as value for the pristine diluting gas; 

(ii) we accept the standard BBN, but assume that lithium 
has been subject to depletion before the cluster formation, 
by means of the reprocessing of the primordial g as in a first 
gener ation of massive, hot stars (population III) l|Piau et al.l 
2006). As in the previous case, we adopt A(Li)=2.20 or 2.3 
in the FG stars and in the diluting FG gas. 

(iii) we assume that Li is depleted in the atmosphere of 
both FG and SG stars, as a consequence of phenomena such 
as diffusion, gravity waves, rotational mixing, or any combi- 
nation of thesfl In this case, we assume A(Li)=2.60 or 2.7 
— as representative of the BBN abundance — for the pristine 
diluting gas. We also assume that the atmospheric effect in 
SG stars is the same that depletes lithium in the FG stars, 
and reduces by <5A(Li)=-0.3 to -0.5 dex the abundance of 
the mixed gas that forms the SG stars. Obviously, a deple- 
tion by 0.3 or 0.5 dex have different weight on the fit of 
observed data. 

In order to begin exploring the range of parameters, we 



Table 3. Possible cases of the lithium dilution when forming 
Second Generation stars and 





Processed ejecta 


Dilution 


Pristine gas 
Abundance 


A 


Supcr-AGB - AGB 


Prist 


BB A(Li)=2.6, 2.7 


B 


Super-AGB - AGB 


Prist 


Pop.II A(Li)=2.2, 2.3 


C 


matter with N(Li)=0 


Prist 


BB A(Li) not relevant 


D 


— A(Li) reduced 


Prist 


like case A 


E 


in Super-AGB — 




like case B 


F 


Supcr-AGB - AGB 


gas lost 








from stars 


N(Li)=0 



2.5 



1 .5 




: .5 



-0.2 



0.2 
[Na/Fe] 



0.4 



0.6 



Figure 2. The lithium versus sodium data in the sample by 
Monaco et al. (2012) for turnoff - subgiant branch stars in M4 
(green open circles with error bars, left scale) are compared with 
the data by D'Orazi & Marino 2010 (scale on the right) for gi- 
ant stars. The abundance in the giants is depleted by convective 
dilution with respect to the abundance in the dwarfs, we assume 
that this dilution is uniform, and compare the data by a simple 
shift of the scale of the A(Li) logarithmic abundances. 



assume an initial A(Li) pr istmc=2.2, or A(Li) pr istmc= 2.7 for 
M4, and A(Li) pr i st inc=2.3, or A(Li) pr istinc=2.6 for the mod- 
els of NGC 2808. This is also very naively motivated by the 
FG abundances listed by different authors for the FG stars 
in different clusters. 



1 We quote an eye-averaged value for the data in the iron range - 
2.5<[Fe/H]<-1.2. Actually, the Asplund et al. (2006) results show 
a clear tren d with metallicity. Th e trend becomes a spread at 
[Fe/H]<-2.7 llSbordone et alj|20ich . 

2 The apparent lack of slope in the Spite plateau is an observa- 
tional constraint that models can achieve by combining diffusion 
with some form of turbulence at the bottom of the atmospheric 
convective zone (Richard et al. 2005; Korn et al. 2006, 2007; Piau 
2008; Lind et al. 2009b). The effect of turbulence is still intro- 
duced in a parametric way in the models. 



2.4 Reference case 

Our standard case studies the SG stellar lithium abun- 
dances predicted by some of the models we presented in 
Paper II and that well reproduce the abundance patterns in 
NGC 2808 and M 4. The lithium abundances in the ejecta 
of super-AGB and massive AGB are taken from Table [2j 
while all the other abundances are listed in Table I of Paper 
II. Table |3] lists as cases A and B the two choices discussed 
above concerning the lithium in the pristine gas. 
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Figure 4. The data for M4 are compared with the simulations, from top to bottom: M4-2, M4-3, M4-4 and M4-5 (Symbols as in 
Fig-0- From left to right, we show 1) case B (A(Li)=2.2 in the pristine diluting gas), 2) case A (A(Li)=2.7); following the two cases in 
which the lithium in the ejecta of the 7.5 and 8 Mq are reduced to A(Li)=2.0, with A(Li)=2.2 (3) or 2.7 (4) in the diluting gas. Finally, 
we show the case in which we assume that the the polluting super- AGB matter has negligible lithium abundance (case C). In this last 
case, the content of the pristine gas (A(Li)=2.2 or A(Li)=2.7) is irrelevant. 



2.5 Other cases 

We compare these "standard" results with two other differ- 
ent possible combinations listed in Table [3] case C in which 
the FG ejecta have no lithium at all (here the initial abun- 
dance of lithium in the pristine gas is irrelevant); cases D & 
E, in which we reduce the yield of the 7.5 and 8 Mq to much 
lower values A(Li)=2.0, and finally we consider a case F, in 
which AGB ejecta are diluted with gas having no lithium. 



3 THE CLUSTERS WITH A MILD O-NA 
ANTICORRELATION 

3.1 M 4 

In Fig . [2] we show the recent data provided by I Monaco et al.1 
(2012) for a sample of main sequence and subgiant stars in 



M 4, and compa re them to the data for g iants in the same 
clusters, given bv lD'Orazi fc Marind (|2010T ). The giants data 
(scale on the right) are shifted to take into account the con- 
vective dilution during the evolution. We notice that the 
two samples provide a slightly different Li-Na trend: the 
M onaco et al.1 (|2012T ) sample shows a slight negative slope 
(the authors show that the difference in A(Li) between the 
low sodium (FG) sa mple and the high sodium ( SG) sample is 
~0.1 dex), while the lD'Orazi fc Marine] |2010l ) data simply 
show a larger scatter in the higher sodium data. The sodium 
zero point also differs by ~0.1 dex between t he two sets. As 
we wi ll compare the models results with the I Monaco et all 
(2012) sample, the slight trend Li-Na is a significant ingre- 
dient of the comparison and affects considerably our inter- 
pretation. 

In Paper I we provided a "standard" model to describe 
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[N,/F»] [N„/F„] 



Figure 3. The sample by Monaco et al. (2011) for turnoff - sub- 
giant branch stars in M4 (open pentagons with error bars) are 
compared with the results of the simulation M4-0 (from Table 2 
by D'Ercole et al. (2012)). Blue triangles are the first generation 
points, black squares are the SG points. The left panel shows case 
B (A(Li)=2.2 in the pristine diluting gas), the right panel repre- 
sents case A (A(Li)=2.7). The red line in both panels represents 
the gas composition along the evolution. The initial part of the 
curve (at large lithium and sodium abundances) is scarcely popu- 
lated due to the fast evolution of the gas composition and to the 
low star formation rate assumed (see Table [TJ. 

the "short" O-Na anticorrelation of clusters like M 4, but 
in Paper II we have shown another class of viable models in 
which SG stars are formed from super-AGB ejecta strongly 
diluted with pristine matter. Fig. [3] shows the results of the 
standard simulation (M4-0) of Paper II, in which the bulk 
of accretion of pristine matter and star formation of the SG 
stars occurs during the evolution of the masses from ~6.5 to 
~5 Mq . In this figure and in the following ones, we display as 
(blue) triangles the FG stars, and as (black) squares the SG 
stars. A stochastic error in the range 0-0.1 dex, compara- 
ble to the observational errors, is introduced, scattering the 
position of the stars, either with respect to the assumed FG 
abundances ([O/Fe]=0.4, [Na/Fe]=0, A(Li)=2.2) for the tri- 
angles, or with respect to the red line, representing the path 
of the gas composition during the SG formation, starting 
from top right and ending to the bottom left. The lithium 
abundance of the SG stars is very similar to that of the FG 
stars, if we are dealing with case B (A(Li)=2.2 or 2.3), be- 
cause the lithium in the ejecta is very similar to what we 
assume to be the abundance in the pristine gas. On the con- 
trary, if we deal with case A (A(Li)=2.6 or 2.7), the abun- 
dances in the ejecta are on average smaller than the pristine 
gas abundance, and the SG shows a slightly reduced lithium 
abundance, in full agreement with the data. 

In Fig. [4] we show instead the results of the simulations 
made in Paper II, in which the super-AGB have a dom- 
inant role for the formation of the SG. It is evident that 
the huge lithium abundances in the 7.5 and 8Mq ejecta 

3 In the standard case for M4 in Paper II we assume that the 
sodium yields for M4 are a factor two larger than those listed in 
Table 1 of that paper. This is done because neon is overabundant 
in the stars of this cluster, so that we should expect that also 
the sodium yield is larger (see the discussion in Paper I). For this 
reason, the two stars with very high lithium in this simulation 
have sodium as large as [Na/Fe]~1.3, a factor two larger than the 
abundance provided by the maximum super-AGB mass 8 Mq 
([Na/Fe]~l). 



have led to an increase in the average SG star lithium abun- 
dance. In the leftmost column of Fig. |4l we show case B 
for cases M4-2, 3, 4 and 5 (from top to bottom). The aver- 
age lithium in the SG is ~0. 15-0.2 dex larger than the FG 
abundance, while the observations show that it is ~0.1 dex 
lower: this full factor two difference rules out case E0. The 
second column from left shows case A: here the situation is 
more ambiguous, and these models can not be excluded (see 
in particular case M4-2). In the third and fourth columns, 
we show respectively case E and case D -where we reduce 
the yield from the 7.5 M and 8 M to A(Li)=2.0: both cases 
are acceptable, apart from the model M4-5 in the last row. 
Finally, the rightmost column shows the results obtained if 
we assume that the ejecta have no lithium, but maintain the 
sodium abundance of AGB-super AGB ejecta in the pollut- 
ing matter (case C). It is obvious that, in this case, assuming 
different Big Bang initial lithium is irrelevant. The results 
show a very good agreement with the observations in all 
cases. 

We see then that the M 4 data do not allow to dis- 
criminate between the AGB scenario and other scenarios in 
which the polluting matter is very sodium rich, but Li-free 
like in the FRMS model. Some models of the AGB scenario 
reproduce well the M 4 Li-Na patterns, but the Li in the 
super-AGB ejecta is not a necessary ingredient. 

Fig. [3]and|4]also show that the models predict the pres- 
ence of a few stars with very large lithium abundance (see 
Sect.[3]for a discussion of the difference in maximum lithium 
between case A and B). In M 4, Monaco et al. (2012) find 
that the star 37934 has A(Li)~2.9 ad relatively high sodium 
([Na/Fe]=0.37). Could it be born in the minute subsample 
of cluster stars formed directly from the super-AGB ejecta, 
just before the strong dilution with pristine gas that gave 
birth to most of the SG? Actually, we notice that the abun- 
dance of sodium predicted by the simulations is too large 
([Na/Fe]~l), but a further exploration of the parameter 
space could be worth. We also remark two consequences of 
this interpretation: 1) this star should be very helium rich. 
Its location at the left border of the main sequence turnoff 
stars (Figure 1 in Monaco et al. 2012) leaves this possibility 
open) ; 2) the HB of M 4 should include a few hot blue tail 
stars, representing the HB phase of these few very helium 
rich stars. 

Another clus ter, NGC 6397 pre sents a case of a super-Li 
rich turnoff star (|Koch et al.ll201ll ). Here A(Li)~4, close to 
our extreme predictions for the 8 Mq super-AGB uncertain 
evolution. According to our models, this super-Li rich star 
should also be characterized by a huge sodium abundance, 
but a preliminary study of its spectrum (private communi- 
cation by K. Lind) appears to indicate a normale sodium 
abundance. 



4 Notice that case B (and, with it, a Big Bang lithium abundance 
A(Li) as low as 2.2), is excluded — if these models represent a 
correct reproduction of the SG formation in M4 — only thanks 
to Monaco et al. (2012) data, showing the mild anticorrelation 
Li-Na. It would not be excluded based on previous observations. 
Thus we remark here how crucial it is to get precise lithium de- 
terminations in multiple populations. 
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3.2 NGC 6397 

The lithium abundances in the stars of NGC 6397 have been 
presented bv lLind et al.l (|2009h . We commented in Paper II 
that NGC 6397 bel ongs to the clusters showing a mild O- 
Na anticorrelation l|Carretta et alj|2009h . and as such is a 
good candidate for a cluster in which the SG is formed by 
super-AGB ejecta strongly diluted with pristine gas. Unfor- 
tunately, we do not yet have models for the super-AGBs of 
metallicity as low as required to discuss more in detail this 
cluster (in particular, we expect smaller oxygen yields), but 
the interpretation of its data can be very similar to that pre- 
sented for M 4 in Fig. [3] A strong dilution of the sodium rich, 
nitrogen rich and oxygen poor super-AGB ejecta may be an 
important factor to explain the problem of th e Be abundance 
in tw o dwarfs of this cluster, analised by iPasquini et al.l 
(2004). Both these stars have beryllium consistent with a 
spallation production of this element, a production that well 
describes the rela tion between beryllium and iron (or oxy- 
gen) in halo stars (|Boesgaard et al.lll999Tl. Nevertheless, one 
of these stars is also oxygen depleted. IPasquini et al.l (|2004r ) 
consider the possibility that this may be a SG star, and cor- 
rectly notice that it should also be beryllium depleted, if 
the polluted matter included in its formation was hot-CNO 
processed. But, if strong dilution of the super-AGB ejecta 
with pristine gas has taken place, the beryllium abundance 
in this SG star may remain close to that of the pristine gas, 
similar to what happens for the lithium abundance. 



4 THE CLUSTERS HARBORING AN 
EXTREME SG 

In this section we focus our attention on clusters that either 
have an extreme O-Na anticorrelation (such as NGC 2808) , 
or for which the presence of a blue main sequence (MS) is ev- 
idence of the presence of a population with very large helium 
content (such as NGC 2808 and u) Cen), or have other signs 
of an high-helium population, like NGC 6752. In this latter 
cluster, the horizontal branch shows a high-temperature tail, 
that puts int o evidence a high-heliu m content in a fraction of 
the SG stars l|D'Antona et al.ll2002T ). a nd the main sequenc e 
too has some indication of splitting (jMonaco et al.l |2010T ) . 
Although this cluster is probably not as extreme as w Cen 
and NGC 2808, it also seems to harbor three MSs identified 
by means of UV colors (Milone et al. 2012, in preparation). 
In addition it ha s an interesting set of O-Li abundances 
jShen et alj|2010h that are worth of an analysis in terms of 
our models. 

4.1 Predictions for the extreme populations: the 
blue MS of NGC 2808 and u Cen 

Paper II reproduces the O-Na anticorrelation and three dif- 
ferent helium contents of the subpopulations in NGC 2808 
by a two-stage process of SG formation, the first episode tak- 
ing place from the pure super-AGB ejecta, leading to the 
very high helium stars populating the blue main sequence, 
the second episode occurring from ejecta diluted with pris- 
tine gas. As discussed, different models are possible for the 
timing and duration of the episode of pristine gas accre- 
tion, all consistent with the observed O-Na anticorrelation. 



In Fig. [5] we now examine the predictions for the Li-Na 
and Li-0 patterns for the model NGC 2808-1 of Paper II, 
that reasonably represents both the O-Na anticorrelation 
and the probable helium difference between the FG and the 
extreme and intermediate SG displayed by NGC 2808. The 
SG points are 40% of the total, following the choice made in 
Paper II. The number of points is proportional to the time 
spent along the line (as well as to the star formation rate) 
so that very few stars have the extreme compositions that 
are linked to the short phase during which the pollutors are 
the most massive super-AGBs. 

4.2 The extreme population 

Large lithium abundances are found in stars formed from 
pure super-AGB ejecta. We notice again the interesting role 
played by the value assumed for the Big Bang (or pris- 
tine) lithium abundance, by comparing the left and right 
panels, representing case B (A(Li) pr j at ine=2.3) and case A 
(A(Li) pr j at ine=2.6). In fact, the maximum Li achieved (for a 
given model of SG formation) is larger in case B, where we 
assume that the pristine gas has the same lithium abundance 
observed in population II stars. In this case, the atmospheric 
abundance coincides with the initial gas abundance, as no re- 
duction is ascribed to mechanisms that deplete it. If the pris- 
tine gas has a larger, e.g. a standard Big Bang A(Li)=2.6, 
we implicitly assume that depletion mechanisms are acting 
at the surface, to bring this value down by ^0.3 dex. These 
same mechanisms will apply also if the lithium abundance 
in the star is the result of the super-AGB ejecta high abun- 
dance. 

Thus, contrary to simple intuition, lithium in SG stars 
(and in particular its maximum abundances, that are di- 
rectly related to star formation from pure ejecta) will be 
~0.3 dex smaller if the diluting gas has a larger abundance, 
in the stars that form from pure ejecta. This could not be no- 
ticed in the bulk of results for M 4 (but is evident in the two 
stars with extreme abundance), due to the strong dilution 
operating in that case. 

An observational confirmation of these predictions 
would provide strong support to the AGB scenario as de- 
scribed by our models. However, we point out that if obser- 
vations will not find a large lithium abundance in the blue 
MS, this would not falsify the AGB scenario but, rather, 
it would simply be an indication that the mass loss rates 
adopted for the models of 7.5 and 8 M©are too extreme (see 
also the discussion of the models for M 4 in section 13. 1[) . 

These observations are not yet available. 

It is important to notice, however, that I Monaco et al.l 

(2010) measured or derived upper limits to the lithium abun- 
dances in 91 MS or early subgiant stars in to Cen, finding 
a remarkable similar value A(Li)=2.19±0.14 for the whole 
sample. In their data, a few stars may have A(Li) close to 
2.4 - 2.5 (still within 3a from the average value), but cer- 
tainly none has a lithium abundance as large as predicted 
in Fig. [5] fo r the blue MS star s in NGC 2808. uj Cen too has 
a blue MS (jBedin et alll2004h th at is interpreted to be very 
helium rich MS (e.g. lNorrisll2004r l. Observations aimed at se- 
lecting specifically blue MS stars to determine their lithium 
abundance could result in some detections of a very large 
lithium content. We remark again that high helium does 
not necessarily mean high lithium, as we show by plotting 
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Figure 5. Li-Na and Li-O simulations for a cluster like 
NGC 2808, harboring an extreme, very helium rich, popula- 
tion and an intermediate one. The parameters refer to model 
NGC 2808-1 in Table 3. We further have A(Li)=2.3 (left panels) 
and 2.6 (right panels). Blue triangles represent the first genera- 
tion population, black squares represent the SG. The red line in 
both panels represents the gas composition along the evolution. 
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Figure 6. Lithium versus helium content in the ejecta of super- 
AGB and AGB stars. The dashed area includes masses of 8, 7.5, 
7, 6.5, 6, 5.5 from top to bottom: all these ejecta would populate a 
"blue MS" , if they are not diluted with helium-poor pristine gas 
before forming SG stars. Only masses from 7.5 to 8 Mq provide 
very lithium rich ejecta. 
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Figure 7. O-Li data for NGC 6752 from Shen et al. 2010. Open 
symbols represent upper limits. In the plot, we have shifted [O/Fe] 
by +0.4 dex, as the data in Shen et al. give a standard [O/Fe]~0 
for the FG of this cluster, contrary to the observations by Carretta 
et al. (2009) that indicate oxygen enhancement for the FG stars 
in this cluster (see text). The line connects the (red) squares that 
show the O-Li yields from Ventura & D'Antona 2009, 2010 and 
2011. 



lithium versus helium content in the ejecta of super-AGB 
and AGB stars in Fig. [6] while masses in the whole dashed 
area (from 5.5 to 8 Mq) may contribute to the blue MS, only 
stars formed from the ejecta of masses in the range 7.5-8 Mq 
would show an outstanding lithium contenlQ 



4.3 The intermediate population 

Another interesting property of the simulations for 
NGC 2808 are the abundances of lithium in the interme- 
diate populations: these result from mixing of pristine gas 
(this latter having a lithium content chosen according to 
cases A or B) with the lithium of the massive AGBs. As 
shown in Fig. [5] the lithium abundance in the intermediate 
population is lower than the pristine abundance, and the 
cases A and B provide interestingly different results. The 
choice of the pristine gas abundance affects the slope of the 
average relation Li-Na and Li-O, that are milder in case 
B. (see the discussion in Sect. 14. 4[) . In any case, the finite 
amount of lithium in the ejecta of massive AGBs produces 
slopes milder than the slope 1 predicted in the case of pure 
dilution of pristine gas with lithium free, sodium-rich and 
oxygen-poor matter. 



5 According to I Monaco et aO 1 120101) , the absence of stars with 
high Li abundance in their oj Cen sample suggests that these were 
not significantly polluted by lithium produced in AGB stars. This 
conclusion is not consistent with the models shown here. 
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4.4 The case of NGC 6752 

NGC 6752 is the first cluster for which the O -Na anticorre- 
lation has been observed in unevolved stars (|Gratton et al.l 
2001), ruling out the possibility that the abundance anoma- 
lies are generated by processes occurring inside observed 
stars, because of the rather low central temperatures and 
thin convective envelopes of stars at the turn-off of GCs. 
Afterwards, this findi ng was confirmed in ot her clusters, 
like M 71 (NGC 683 8 iRami'rez fc Cohenll2002h and 47 Tuc 
|Carretta et~ai1l2004l ) Further data for gian ts of NGC 6752 
have been collected in lCarretta et all (|2007l ). The data show 
a typical O-Na anticorrelation, not extremely extended in 
oxygen, although several upper limits to [O/Fe] are present. 
The sodium-poor stars (FG stars) have largely a-enhanced 
[O/Fe]~0.4-0.5. 

Lithiu m in NGC 6 7 52 w as measured in nine turnoff 
stars by IPasquini et alj (|2005f ). who found that its abun- 
dance was anticorrelated with sodium and nitrogen, and di- 
rectly correlated with oxygen. These data formed the basis 
of an interpretation based on a dilution model of pristine 
matter with the matter processed i n massive rotating stars , 
having n o lithium an d high sodium l|Decressin et al.ll2007bl ) . 
Recently. IShen et al.l (|2010h examined a much larger sample 
(112 stars), confirming the O-Li correlation, but finding a 
slope of the correlation AA(Li)/A[O/Fe]~0.4, much milder 
than the slope 1 expected if the lithium in the pristine mat- 
ter were just diluted with Na-rich Li-free ejecta. Figure [7] 
shows the O-Li data for NGC 6752 together with our yields 
for super-AGB and AGB stars. Two caveats must be no- 
ticed: 1) here we plot the computed abundances, as we deal 
with turnoff stars, in which no kind of further deep mixing 
to reduce oxygen can be invoked; 2) in the comparis on with 
simulations, we have shifted the en et all (|2010f ) [O/Fe] 
abundan ces by +0.4 dex: in fact, the "standard" [O/Fe] 
values in IShen et al l |201Gft analysis is ~0. 0-0.1, implying 
some systemat i c shif t with respect to the abundances of 
ICarretta et all l|2007t) . Our initial models have [O/Fe]=0.4, 
so the observational values have been adjusted. 

The O-Na data for this cluster are not as good as the 
data for NGC 2808, and also the information of the MS 
are still incomplete in the literature. Therefore, we use as a 
guideline the simulations made for NGC 2808. In Fi g. [5] we 
overplot the simulations of Fig. [5] to the IShen et al.l (|2010l ) 
data, after shifting the data by +0.4dex in [O/Fe]. We also 
shift A(Li) of the FG by +0.05 dex to achieve a better re- 
production of the undepleted abundances in NGC 6752. We 
adopt the same shift in oxygen for the SG simulation, as the 
oxygen depletion is directly dependent on the initial abun- 
dance, but we do not shift the SG lithium upwards, because 
the yields in lithium are not linked to the initial abundance 
chosen for the models, but only to the lithium production 
through the Cameron Fowler mechanism during the HBB. 
Fig. [S] shows that the abundances in the ejecta may play a 
role in decreasing the slope of the [0/Fe]-A(Li) anticorrela- 
tion. There is no sign of the few very Li-rich stars predicted 
for the extreme second generation of NGC 2808, but we 
see that the "intermediate" population is broadly consistent 
with the data. By comparing the two panels we see that 
the slope of the simulated data points (excluding the super- 
AGB values) differs in the two cases (A & B): as expected, 
the slope is flatter in the case of A(Li)=2.3. We compute an 



average slope of the simulation by a least square fit both of 
the FG and SG points, after excluding the extreme lithium 
abundances. We obtain AA(Li)/A[O/Fe]=0.27±0.05 for the 
case A(Li)=2.3 and AA(Li)/A[O/Fe]=0.62±0.05 for the 
case A(Li)=2.6. Thus further observations of the lithium 
patterns in FG and SG stars providing a more robust de- 
termination of the [0/Fe]-A(Li) slope, along with a more 
detailed model in which also other chemical constraints are 
considered, may even result in a test of the standard BBN. 



5 MODELS IN WHICH THE DILUTING GAS 
DOES NOT CONTAIN LITHIUM 

We finally shortly examine case F (Table 3), in which the di- 
luting matter is devoid of lithium. In order to solve the prob- 
lem of the source of the matter required for dilution, Grat- 
ton & Carretta (2010) proposed that it comes from mass 
loss from FG stars, that can quantitatively provide the gas 
necessary to achieve this, if we consider the main sequence 
stars of all masses and reasonably large mass loss rates. Con- 
sequently, this gas has no lithium, because lithium burns be- 
low the stellar surface as soon as the temperature increases 
above ~ 2 x 10 6 K, and the mass loss rate required in this 
scenario implies that layers of matter that have reached this 
temperature must contribute to the dilution. Another pos- 
sible case is when the diluting gas comes from the mass 
lost by non conservative evolution of massive close binaries 
jVanbeveren et all I2OI2T ) or by close encounters of stars in 
the dense cores of the GCs (Carini et al. 2012). It is clear 
that the models for M4 in which there is a strong dilution 
with pristine matter (from M4-1 to M4-5) lithium in the SG 
will be very small, if any, in contrast with the observations. 
As for NGC 2808, in Fig. [9] we show the implications of a 
lithium-free diluting gas on the Li-O and Li-Na patterns in 
the model 2808-1. Contrary to the cases shown in Fig. 5, this 
simulation predicts much smaller lithium abundances for the 
intermediate population. Further observations of lithium in 
clusters may help to falsify this model. 



6 THE BIG BANG LITHIUM ABUNDANCE 

This work, first aimed at further discussing a detailed chemi- 
cal evolution model for the multiple populations of GCs, has 
revealed an interesting implication for the Big Bang nucle- 
osynthesis. In fact, despite the scarcity of lithium abundance 
determinations in GCs, these data may result to be a power- 
ful and independent way to constrain the lithium abundance 
emerging from the Big Bang. The abundance of lithium in 
the gas forming the SG stars contains information on the 
abundance of lithium in the gas in which the hot-CNO pro- 
cessed ejecta are diluted. If we can assume (and this is very 
reasonable) that, at first order, the possible phenomena that 
deplete lithium at the surface of population II stars are in- 
dependent of their being FG or SG stars, the difference in 
lithium between FG and SG stars will keep track of the pris- 
tine gas abundance. 

In the case of M4, we have seen that the cases B 
(A(Li)=2.2 in the pristine gas) shown in Fig. 4 can be ex- 
cluded. This could be a very interesting result, but it is 
subject to a confirmation of the large lithium yields shown 
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Figure 8. Overplotted to the Li-O for NGC 6752 (shifted by +0.4dex, as in Fig.O, we show the results of simulations of Fig. [5] Green 
squares represent the FG and green triangles the SG. The initial lithium abundance of the diluting pristine gas is A(Li)=2.6 (left panel) 
and 2.3 (right panel). To achieve a better fit of the FG data, we have shifted A(Li) of the FG by +0.05 dex. 
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Figure 9. Li-Na and Li-O simulations for model NGC 2808-1 in 
Table 3. The difference with respect to Fig. 5 is in the assumption 
of no lithium in the diluting gas. We further have A(Li)=2.3 (left 
panels) and 2.6 (right panels). Notice that the extreme population 
is similar to that shown in Fig. 5, because the SG formation occurs 
in the undiluted matter, while the intermediate population, where 
dilution is important, has much lower lithium abundances. 

in Table 2 for super-AGB stars. As far as we have seen from 
the other observations, there is no present confirmation of 
those yields. Further, they are theoretically very uncertain, 
being linked to the choice (unconstrained by observations) of 



the mass loss rate adopted in these computations. A confir- 
mation could come from observation of large lithium abun- 
dances in the blue MS stars of NGC 2808 or uj Cen. If in 
the end the yields of super-AGB stars will be constrained 
to be much smaller than those used here, case B will still be 
possible, and we will have to resort to different comparisons. 

A very interesting information seems to come from the 
slope of the O-Li data in NGC 6752: Figure 8 shows the 
possibility of constraining the primordial lithium abundance 
from the data. This result is less dependent on the super- 
AGB yields, as it refers to the (later) phase of formation 
of the intermediate population. Of course, more and better 
data, and stronger constraints on the dilution model and 
on the AGB yields are necessary. We only wish to point out 
this rewarding possibility, to encourage further and extensive 
observations of lithium abundances in clusters with multiple 
populations. 



7 DISCUSSION 

In this paper we have presented the first predictions for the 
abundance of Lithium in first and second generation stars of 
GCs, based on quantitative models that can succesfully re- 
produce the O-Na anticorrelation and the helium abundance 
distribution in the clusters M 4 and NGC 2808 (Paper II). 
Although the model is necessarily parametric, it allows us 
to discuss the lithium data appeared in the recent literature 
for several GCs: NGC 6752, M 4, NGC 6397, u Cen. Our 
models are based on our recent yield predictions for massive 
AGB and super-AGB stars. We show that: 

(i) An important ingredient in the study is the dilution of 
the ejecta with the pristine gas, and thus whether we inter- 
pret the surface abundance of lithium in population II stars 
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as the primordial abundance, or as due to in situ depletion 
mechanisms. A better understanding of the SG formation 
could perhaps help to constrain better the Big Bang abun- 
dance, although this is at present only wishful thinking. 

(ii) Not all models for M 4 are consistent with the slight 
decline in lithium by 0.1 dex between the FG and SG. A 
good fit of the data is achieved also assuming zero lithium in 
the super-AGB matter (rightmost column panels in Fig[4]), 
thanks to the large sodium yields of these ejecta, that require 
strong dilution with pristine matter. Models in which the 
O-Na "short" anticorrelation is the result of the rapid for- 
mation of SG stars from strongly diluted super-AGB ejecta 
wo uld solve the mystery of the Be-rich, O-poor star found 
bv lPasauini et all l|2004 ) in NGC 6397. 

(iii) In complex clusters, in which a very He-rich extreme 
population is present, we should expect that this population 
is born from pure super-AGB ejecta, and thus is very Li- 
rich. Observations of lithium in the blue MS of w Cen and 
NGC 2808 may verify or confute this prediction. 

(iv) The O-Li observations of NGC 6752 are consistent 
with the intermediate population of the models discussed 
for NGC 2808. The slope of the correlation (milder than 
the slope 1 expected if the polluting matter is Li-free) is 
reproduced, thanks to the finite lithium content of the AGB 
ejecta. In this case too, a better determination of this slope 
might provide hints to the Big Bang Lithium abundance. 

(v) In some models for the SG formation, the diluting 
gas comes from nuclearly unprocessed matter from FG stars, 
that is probably Li-free. We predict that in this case lithium 
in SG stars should be much smaller than in the FG. 
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